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By combining experimental measurements of the quasiparticle and dynamical magnetic proper- 
ties of optimally electron-doped Pro.88LaCeo.i2Cu04 with theoretical calculations we demonstrate 
that the conventional fermiology approach cannot possibly account for the magnetic fluctuations 
in these materials. In particular, we perform tunneling experiments on the very same sample for 
which a dynamical magnetic resonance has been reported recently and use photoemission data by 
others on a similar sample to characterize the fermionic quasiparticle excitations in great detail. 
We subsequently use this information to calculate the magnetic response within the conventional 
fermiology framework as applied in a large body of work for the hole-doped superconductors to find 
a profound disagreement between the theoretical expectations and the measurements: this approach 
predicts a step-like feature rather than a sharp resonance peak, it underestimates the intensity of 
the resonance by an order of magnitude, it suggests an unreasonable temperature dependence of the 
resonance, and most severely, it predicts that most of the spectral weight resides in incommensu- 
rate wings which are a key feature of the hole-doped cuprates but have never been observed in the 
electron-doped counterparts. Our findings strongly suggest that the magnetic fluctuations refiect 
the quantum-mechanical competition between antiferromagnetic and superconducting orders. 

PACS numbers: 74.72.-h, 74.25.Ha, 74.20.Rp, 75.40. Gb 



I. INTRODUCTION 

High-Tc superconductivity occurs in doped Mott- 
insulatorsi, electronic states which are insulating because 
of dominating electron-electron interactions and charac- 
terized by spin-only antiferromagnetism. The study of 
the fate of this magnetism in the doped, superconduct- 
ing systems has been on the forefront of high- research 
from the very beginning with inelastic neutron scatter- 
ing being the primary source of experimental informa- 
tion. This is now well documented in the hole doped 
(p-type) superconductors where one finds the famous 
" hourglass" spectrum of magnetic fluctuations^i^i^i^i^iii^, 
but the interpretation of these findings has been subject 
of severe controversy. On the one hand this can be in- 
terpreted as the signature of strong interaction physics, 
where the magnetic fluctuations signal the competition 
between superconductivity and incommensurate, Mott- 
like antiferromagnetism (the " stripes" what mat- 

ters most in this interpretation is that it is envisaged 
that the measurements reflect the quantum dynamics 
of competing, highly collective order parameter fields^^. 
However, it turns out that these data find also a credi- 
ble interpretation in terms of conventional Fermi-liquid 
(FL) physicsi^ii^ii^ii^. Here, it is asserted that the elec- 
tron system renormalizes in a weakly interacting Fermi- 
gas, acquiring a conventional Bardeen-Cooper-Schrieffer 
(BCS) gap in the superconducting state. The residual 
interactions then give rise to weakly bound states resid- 
ing in the gap formed in the non-interacting particle-hole 



spectrum, governed by the random phase approximation 
(RPA). 

Only very recently inelastic neutron scattering data be- 
came available for the magnetic fluctuations in the elec- 
tron doped (n-type) superconductor Pro.88LaCeo.i2Cu04 
(PLCCO)^-. The spectrum is dominated by a dynamical 
peak (resonance) at an energy lOj-cs ~ II meV residing at 
the antiferromagnetic wavevector qaf = {t^, tt) whereas 
the incommensurate branches (wings) found in the p- 
type superconductors in the vicinity of the resonance are 
conspicuously absent. 

Here we will employ tunneling spectra obtained for 
the same sample as used for the neutron measurement, 
in combination with angular resolved photo emission 
spectroscopy (ARPES) by others on a similar sample 
to characterize the fermionic quasiparticle excitations 
in great detail. We subsequently use this information 
to derive the magnetic spectrum employing the RPA, 
to find out that there is a profound disagreement be- 
tween the theoretical predictions for the magnetic fluctu- 
ations coming from this fermiology interpretation and the 
measurementslS.. In particular, (i) this framework pre- 
dicts a very asymmetric almost step-like feature slightly 
above the edge of the particle-hole continuum instead of 
a sharp resonance peak seen in neutron scattering, (ii) 
it suggests a strong temperature dependence of the res- 
onance feature, both in intensity and position, inconsis- 
tent with the data, (iii) it underestimates the absolute 
intensity of the resonance by an order of magnitude, and 
finally (iv) it predicts that most of the spectral weight re- 
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sides in incommensurate wings below the resonance fea- 
ture in clear contradiction to the data^^. 

The outline of this paper is as follows. In Section HTl we 
explain the workings of the FL/RPA approach and the 
extraction of the quasiparticle parameters from ARPES 
data and our tunneling experiments. The latter we de- 
scribe in detail in Section IIIII The results of the theo- 
retical calculations are presented in Section HVl and com- 
pared to the magnetic excitation spectrum of PLCCO. 
Finally, our results and implications of our findings are 
summarized and discussed in Section fVl 



Xo(q, 1^), respectively. The latter condition, the vanish- 
ing of the bare Lindhard function Xo(q, w), enforces a 
resonance at (q, uj) to live at an energy uj below the gap 
of the particle-hole continuum, lo < 2A(q). 

Before we can employ the RPA formula ([T]) to calculate 
the magnetic response x(qj'^)j we have to characterize 
the bare quasiparticles in great detail to determine the 
bare BCS susceptibility In particular, we have to use 
experimental input to extract the normal state dispersion 
e(q) and the superconducting gap A(q). For the normal 
state dispersion we use the standard square lattice tight- 
binding dispersion 



II. DETAILS OF THE FL/RPA CALCULATIONS 

Let us first describe the standard calculations based on 
the FL/RPA framework. In this approach it is assumed 
that the cuprates can be interpreted as FLs all along 
(including the normal state) undergoing a weak coupling 
BCS instability towards a d-wave superconductor, while 
the excitations are calculated from the leading order in 
perturbation theory (RPA) controlled by the weakness of 
the residual interactions. The spin susceptibility within 
RPA can be written as 



x(q,^) 



xo(q,t^) 



1 - U{q)xo{n,t^) 



(1) 



where U{q) denotes the fermionic four-point vertex and 
Xo(q, t^) the bare non-interacting BCS susceptibility, 
which is completely determined by the normal state 
tight-binding dispersion e(q) and the superconducting 
gap function A(q), namelj^ 
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Here £'(q) = Y^e^(q) + A^(q) denotes the quasiparticle 
dispersion, / the Fermi function, and for abbreviation 
we have defined fik.q = (ek-i-qCk + Ak+qAk)/(-Bk+q-Bk)- 
The three parts in xoi^j^) are due to quasiparticle scat- 
tering, quasiparticle pair creation and quasiparticle pair 
annihilation, respectively. 

In the FL/RPA approach for the magnetic resonance 
mode of the p-type cuprates, the dispersing incommen- 
surate wings merging into the commensurate resonance 
peak at qAp are interpreted as a dispersing bound state 
formed in the gap below the particle-hole continuum. 
Such a bound-state corresponds to a pole in the imag- 
inary part of the susceptibility, x"(q, w), given by the 
conditions 1 — U{q)x'o{'i, w) = and Xo (q^ t^) = for the 
real and imaginary part of the bare BCS susceptibility 



e(q) = —2t[cos{kx) + cos{ky)] ~ 4:t' cos{kx) cos{ky) 
-2t"[cos{2kx) + cos(2fcj^)] 
— 4i"'[cos(2fcj,) cos(fcj^) + cos{kx) cos(2fcy)] 
-4f " cos(2fc^) cos(2fcj,) - fj,, (3) 

having incorporated an appropriate chemical potential /i. 
A normal state single particle dispersion e(q) for opti- 
mally doped PLCCO of this form has been determinedi^ 
by fitting the ARPES data at 30 along three indepen- 
dent directions. The resulting tight-binding parameters 
are listed in Table [D 
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TABLE I: Collection of parameters used in our calcula- 
tion: parameters of the normal state tight-binding disper- 
sion t,t' ,t" ,t"' ,t^^ , chemical potential /i, d-wave gap parame- 
ters Ai, A3, and four point vertex parameters U,AU. Tight- 
binding parameters for PLCCO are taken from Ref. |3 and 
parameters for YBCO from Ref. [T^ . 



To reproduce the non-monotonic d-wave gap of 
PLCCO observed in the ARPES measurement's,, we in- 
clude third harmonics in the gap function, 



A(k) 



^[cos{k^) - cos(/cy)] 



A3 
2 



C0s(3fca;) — C0S(3fc.y)] 



(4) 



and adjust the ratio A1/A3 to reproduce the functional 
form of the gap along the Fermi surface found experimen- 
tally. For A1/A3 « 2.43 we find a maximum gap value 
Amax ~ l-3Ao under a Fermi surface angle 0max ~ 21° 
with Aq the gap value at the anti-nodal direction (0 = 0) 
in agreement with the experimental observation (see Fig. 
[1]). The gap maxima are very close to the intersection 
points of the Fermi surface and the magnetic Brillouin 
zone I fca; I -I- 1 fcj, I < TT. These so called hot spots are relevant 
for particle-hole processes contributing to the magnetic 
response at qAp- 
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FIG. 1: The non-monotonic d-wave gap A of PLCCO along 
the Fermi surface as a function of the Fermi surface angle (p 
(see inset) calculated with the set of parameters listed in Table 
IT] The inset shows the relation between the Fermi surface and 
the magnetic Brillouin zone. The hot spots relevant for the 
magnetic response at qaf are shown as open circles. The po- 
sition of the gap maximum close to the hot spots and the ratio 
Amax/Ao of the maximum gap value and the antinodal gap 
are in good agreement with ARPES measurements^^. The 
absolute gap values are extracted from our tunneling experi- 
ment (see Section . 



To determine absolute gap values which are difficult 
to extract from leading-edge shifts in ARPES data we 
have performed tunneling experiments on the same sam- 
ple of PLCCO showing a magnetic resonance at uj^cs ~ 11 
meV— . A detailed discussion of the experimental setup 
and the obtained results are presented in section Hill 

The complete set of parameters for the normal-state 
dispersion e(q) and the gap function A(q) of PLCCO is 
listed in TableUand compared to a set of parameters used 
recently to calculate the magnetic response of optimally 
doped YBCOi^. The latter we use as a benchmark for 
our numerical calculation and also for a comparison of 
the features of the FL/RPA spectra of the n- and p-type 
compounds. 

To calculate the bare susceptibility xoi^,^), we re- 
place iO^ by iF in the energy denominators, mimicking 
experimental broadening. We take F = 2 meV consis- 
tent with the typical broadening in neutron scattering 
and values used in other RPA calculations. The result- 
ing well-behaved function is then summed numerically 
over a 1500 by 1500 mesh in the Brillouin zone. 

Since the bare non-interacting BCS susceptibility 
Xoi^.,'-^) is completely determined by e(q) and A(q), we 
can only adjust U (q) in the RPA equation ([T]) to repro- 
duce the magnetic excitation spectrum of PLCCO. Fol- 
lowing other standard RPA calculations for p-type com- 
pounds, we take an onsite repulsion Uq and allow for a 
small q-modulation with amplitude AU (see e.g. Ref. 
M), f/(q) = Uo- AU[cos{q,) + cos{qy)]. 



III. TUNNELING EXPERIMENT 



To determine the absolute gap value of PLCCO and its 
temperature dependence, we performed tunneling mea- 
surements on the same sample used for the neutron scat- 
tering measurements^!. The directional point-contact 
tunneling measurements were carried out by pointing a 
Au tip towards the specified directions of a or b crystal 
axis which is determined by neutron scattering (Fig. [2] 
a,b). The Au tips were mechanically sharpened by care- 
fully clipping a gold wire with a diameter of 0.25mm. 
The approaching of the tips were controlled by a refined 
differential screw. The point contact insert was set in the 
sample chamber of an Oxford cryogenic system Maglab- 
EXA-12. In order to reduce the quasiparticle scattering 
in the barrier layer and hence obtain high quality data, 
the nonaqueous chemical etch was used to attenuate the 
insulating layer on the sample surface immediately before 
mounting the sample on the point contact device^S. 

Typical four-terminal and lock-in techniques were used 
to measure the I V curves and the differential resis- 
tance dV/ dl vs V of the point contacts simultaneously. 
Then the dynamical conductance dl/dV ~ V was ob- 
tained both by converting the dV/ dl ^ V curves and by 
calculating the derivative of / ~ ^ relations in order to 
ensure the reliability of the results. It was verified that 
the results were not affected by the heat-relaxation effect 
by comparing the curves recorded by positively and nega- 
tively bias scanning. For quantitative analysis, the spec- 
tra were normalized by corresponding backgrounds con- 
structed according to the spectrum measured well above 
T 

In Fig. [2ji, we show the temperature dependence of 
the dl/dV spectra from 2 K to 20 K with increments of 
2 K. Note that due to the experimental broadening the 
two van-Hove singularities at Aq and A^ax (see Fig|T]) in 
the density of states p{uj) (Fig. 2c) are not resolved. To 
make it as advantageous as possible for the the FL/RPA 
approach to explain the magnetic resonance, we iden- 
tify the gap seen in the tunneling spectra with the gap 
Aq at the antinodal direction. This probably overesti- 
mates the true gap since from the data we probably ex- 
tract an energy between Aq and A^ax- On the other 
hand, we note that point-contact tunneling measures the 
density-of-states averaged superconducting gap, its value 
might be different from those obtained by spatially re- 
solved scanning tunneling microscopy. 

From a fit to the extended Blonder-Tinkham-Klapwijk 
(BTK) modeP^ with a d-wave- type gap function^", we 
obtain the BCS like temperature dependence of the gap 
value as shown in Fig. [2^. Similarly, from the depen- 
dence of the spectra on c-axis aligned magnetic field, we 
extract the superconducting gap as a function of increas- 
ing magnetic field (Fig. [lf,g). 
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FIG. 2: Geometry and results of direct point-contact tun- 
neling measurements on single crystals of PLCCO. (a) The 
schematic diagram of the experimental setup, where a An 
tip is pointed towards the a/b axis direction determined by 
neutron diffraction, (b) The relationship between standard 
d-wave gap and tunneling direction, (c) Calculated quasipar- 
ticle density of states using gap values at different tempera- 
tures showing Van Hove singularities at the antinodal gap Ao 
and the maximum gap Amax. (d) Temperature dependence of 
the dl /dV spectra from 2 K to 20 K every 2 K. The spectra 
were obtained by normalizing the corresponding backgrounds 
at temperatures well above Tc- (e) Temperature dependence 
of the gap value Ao, the solid line denotes the BCS predic- 
tion, (f) Magnetic field dependence of the dl /dV spectra for 
a c-sxis aligned magnetic field. The theoretical calculations 
are indicated by red lines in (d) and (f ) , respectively. All the 
spectra and fitting lines except for the lowest ones are shifted 
upwards for clarity, (g) Superconducting gap as a function of 
increasing magnetic field, the solid line is the guide to eyes. 
The Ao values in (e) and (g) are determined by fitting the nor- 
malized spectra to the extended Blonder- Tinkham-Klapwijk 
model^i with a d- wave- type gap function along the a/h axes. 



IV. RPA RESULTS AND COMPARISON TO 
EXPERIMENTS 

Before we calculate the magnetic response x"(<Ij ^) fo^' 
PLCCO within the FL/RPA framework using the tight- 
binding dispersion e(q) and the gap function A(q) deter- 
mined by ARPES— and our tunneling experiment, we 
first test our numerical routine for a set of parameters 
that has been used to calculate the magnetic excitation 



spectrum of optimally doped YBCOli. The resulting 
magnetic excitation spectrum in the vicinity of the anti- 
ferromagnetic wave vector is shown in Fig. [3] along the 
{H, 1/2) and [H, H) directions, respectively and is found 
to be ill perfect agreement with the theoretical results 
in Ref. [Tj. The favorable comparison of theoretical re- 
sults with the dispersion found in inelastic neutron scat- 
tering experiments'^'^ on optimally doped YBCO is also 
shown in Fig. [3l However, a closer inspection of the 
intensities shows that the FL/RPA calculation severely 
underestimates the spectral weight above the commen- 
surate dynamical resonance^^. Whereas experimentally 
the intensities of the upper and lower wings forming the 
characteristic hourglass in the vicinity of the resonance 
are quite comparable^ in the RPA results the upper half 
of the hourglass is completely absent (see Fig. . 



q=(H,H) (r.l.u.) 

0.3 0.5 0.7 




0.3 0.5 0.7 

q={H,0.5) (r.l.u.) 



X"(q.ra) (meV)-i 

FIG. 3: Magnetic response calculated within the 

FL/RPA approach using a set of parameters optimized for op- 
timally doped YBCOi^ plotted along the {H, 0.5) and {H, H) 
(inset) directions close to qaf. White and black points show 
neutron scattering data from Ref. and Ref. d, respectively. 



A. Resonance feature of PLCCO 

Before we calculate the full momentum dependent 
RPA spectrum of PLCCO, we try to reproduce the reso- 
nance feature at io^cs 

w llmeV found in inelastic neutron 
scattering^^ by tuning the value of the four-point vertex 
U — U{(\af) at the antiferromagnetic wavevector. 

The smallness of the gap A « 5meV relevant for qaf- 
scattering (see Fig. [T]) enforces a worrisome fine-tuning 
to produce a bound state. The necessary conditions for 
the corresponding singularity in the imaginary part of the 
dynamic susceptibility x"(qAF,^) are given by w < 2 A 
and U = l/xo(qAF, ^^)- The evolution of x"{^) for dif- 
ferent values of U is plotted in Fig. [H For U < 515meV 
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the resonance is pushed into the particle-hole contin- 
uum whereas the system runs into a magnetic instabil- 
ity for U > 528meV. Since experimentally the resonance 
peak is found slightly above the edge of the particle-hole 
continuum as confirmed by our tunneling measurements 
(Fig. Efcjg) it cannot be explained as a bound state. For 
U = 500meV, we find an intensity enhancement around 
llmeV. As expected, since no bound state is formed the 
FL/RPA result has a very asymmetric and almost step- 
like line-shape, rather than the symmetric peak observed 
in experiment^ (see Fig. |6b), and the intensity is signif- 
icantly reduced compared to a typical bound state situ- 
ation. 
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FIG. 4: Evaluation of x" {<Iaf,ui) in the superconducting state 
for different values of U showing the narrow energy window 
515meV< U < 528meV for which a bound state in the gap 
2A — lOmeV of the particle-hole is formed. For U — 500meV 
we find an intensity enhancement at the experimentally ob- 
served resonance energy Wres = llmeV. Since the feature is 
located slightly above the gap the intensity is significantly 
reduced compared to the bound-state situation and the line- 
shape is very asymmetric. 



B. Temperature dependence 

In this section we are going to analyze the tempera- 
ture dependence of the resonance feature that is to be 
expected in the FL/RPA framework taking the BCS like 
temperature dependence of the gap given by the tunnel- 
ing experiment. Since the two features in the quasipar- 
ticle density of states at the antinodal gap Aq and the 
maximum gap Amax are not resolved in the data (see 
Fig. ^,d) we assume the non-monotonic functional form 
of the gap along the Fermi surface not to change with 
temperature and simply scale the gap function plotted 
in Fig. [1] according to BCS like temperature dependence 
extracted from the tunneling data (Fig. 

The resulting temperature dependence of the "res- 
onance" feature we obtained at zero temperature at 
llmeV for U = 500meV is summarized in Fig. [S] From 



the calculation, below Tc « 24K we expect a strong 
temperature dependence of the resonance feature both 
in position and intensity. With increasing temperature 
the resonance shifts to lower energies whereas the inten- 
sity goes down continuously (Fig. [5^,b). These predic- 
tions are inconsistent with the experimental observations, 
where the position of the resonance appears to be fixed 
and the intensity drops down sharply close to Tj^. 
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FIG. 5: (a) Evolution of the resonance feature calculated 
for PLCCO within the FL/RPA approach with temperature. 
With increasing temperature the resonance feature shifts con- 
tinuously to lower energy and decreases in intensity, (b) Cal- 
culated intensity at a function of temperature showing 
a strong decrease at temperatures well below Tc. 

Since the temperature dependence of the gap is ex- 
pected to be the dominant effect we have not taken ther- 
mal broadening into account. However, additional broad- 
ening would even lead to a stronger continuous decrease 
of the intensity below Tc- 



C. Comparison of absolute intensities 

Since the magnetic resonance of PLCCO is too high in 
energy it cannot be explained as a bound state within the 
FL/RPA framework and has consequently to be identi- 
fied with a weak intensity enhancement slightly above the 
gap of the particle-hole continuum. Therefore we expect 
a significant reduction of the spectral weight compared 
to a typical bound-state situation. Comparing the inten- 
sity of the "resonance" feature we produced within the 
RPA calculation for U = 500meV at Wrcs = llmeV with 
the intensity of the bound state one obtains with the set 
of parameters optimized for YBCO^^ (see Table we 
expect the resonance of PLCCO to be weaker by a factor 
15 than the resonance of YBCO. In this comparison we 
have used the same broadening F — 2meV for both cases. 

To compare this theoretical expectation with experi- 
ment, we have converted the neutron scattering raw data 
on the resonance of PLCCO reported in Ref. \1Q to abso- 
lute units (^^eV~^f.u.~^) both in the normal and su- 
perconducting states by normalizing them to acoustic 
phonons around the (2,0,0) Bragg refiectioni. 

In the long- wavelength limit, the differential cross sec- 
tion for coherent one phonon emission at given (/?, uj) is^ 
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where k = t + q is the momentum transfer of the neu- 
tron, E{q) the energy of the phonon mode, N the num- 
ber of unit cells, ki and kf are the incident and final 
wavelengths of the neutron, n(w) is the standard Bose 
population factor, Cqs is the unit vector in the direction 
of atomic displacement for the phonon mode, e~^^ is 
the Debye- Waller factor, M the mass of the unit cell, 
and G(t) is the standard nuclear structure factor. The 
spectrometer dependent constant A can be determined 
through the measurement of a known phonon in the ma- 
terial. For our case, we measured a transverse acoustic 
phonon at Q = (0.12,2,0). 

The same spectrometer dependent constant A can then 
be used to determine the magnetic susceptibility in abso- 
lute units. For paramagnetic spin fluctuations the cross 
section is 



dndE 
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where (7ro)^/4 is 7.265 ■ 10~^^ cm^ and /(k) is the 
isotropic, magnetic form factor for Cu^^. In order to ob- 
tain the local susceptibihty x"i^) — J x"(Qi^)d"^Q 
at the (tt, tt) in-plane wavevector, Q-scans were per- 
formed at selected energies. For energies below 5 meV 
and above 10 meV, SPINS data and BT-9 data were 
respectively cross-normalized to the absolute values of 
the HB-1 data using constant scale factors. For ener- 
gies below 5 meV, the measured Q- widths along [H, H] 
were broader than resolution while scans at all higher en- 
ergy transfers all showed resolution limited peaks along 
[H, H] . In order to estimate the local susceptibility, the 
magnetic signal was assumed to be a two-dimensional 
Gaussian within the [H, K] plane and rod-like out of 
plane. This neglects the rotation of the resolution el- 
lipsoid at energy transfers away from the resonance po- 
sition and results in a slight underestimation of the in- 
tegrated magnetic scattering at energies below the reso- 
nance. This estimation however is systematic and does 
not influence relative changes in the local susceptibility as 
the system enters the superconducting phase. For points 
in i?-scans with E > 5 meV in which no Q-scan data was 
available, the calculated resolution value was used pro- 
jected along the [H, H] direction. The background was 
removed through subtracting the measured nonmagnetic 
signal away from the correlated (tt, tt) position as shown 
in Ref. [3. All data was corrected for A/2 contamination 
in the monitor, and in our calculations for data at both 
2 K and 30 K, the Debye- Waller factor was assumed to 
be 1. 



Assuming all the scattering centered at Q = 
(1/2,1/2,0) is magnetic, we find that the local suscep- 
tibility x"{lli) has a peak around 11 meV and increases 
at all energies probed (from 0.5 meV to 16 mcV) on cool- 
ing from the normal state to the superconducting state 
(Fig. EJd). This is in contradiction to the theoretical 
results which predict a reshuffling of spectral from low 
energies in the normal state to the resonance feature in 
the superconducting state (Figs. [5^(6)3). 

Figure [B^ shows the local susceptibility for optimally 
doped YBCO'*, where the resonance intensity is obtained 
by taking the temperature difference between the normal 
(100 K) and superconducting states (10 K) since the ab- 
solute intensity of the mode in the normal state is still 
unknown. In Figure [BJa we plot the local susceptibility 
in the normal (30 K) and superconducting (2 K) states 
normalized to phonons. The local susceptibility in abso- 
lute unit is similar to those of PLCCO with a different 
Tc^^ and Pro.89LaCeo.iiCu04^^, and is about 2.5 times 
smaller than that of the resonance for YBCO in Figure 
|6^. From the FL/RPA calculations we expect the spec- 
tral weight of the resonance to be smaller by a factor 15 
as compared to YBCO's resonance (see Fig. [S]), in clear 
contrast to experiments. 



D. Momentum dependence 

Finally, we calculate the momentum dependence of the 
imaginary part x"(q, w) of the dynamic susceptibility in 
the vicinity of qaf using the band structure parameters 
and superconducting gap discussed in Sec. |lT] as appro- 
priate for PLCCO. In Sec. IIV Al we have seen that for 
[/(qaf) = 500 meV the FL/RPA approach reproduces a 
feature at Wres = 11 meV. However, since this feature is 
located at an energy above the gap of the particle-hole 
continuum its lineshape and spectral weight are inconsis- 
tent with the experimental observations. 

We start with a momentum independent four-point 
vertex U{q) — U (Hubbard- like approximation) which 
in the the case of p-type compounds turns out to give 
a pretty good description of the magnetic excitation 
spectra^ii^, much better than for a strong momentum 
dependent J7(q)i^. 

Using a constant U — 500 meV producing a resonance 
feature at the experimentally observed energy Wres = 11 
meV^- the FL/RPA predicts for the n-type supercon- 
ductor the spectrum shown in Figure [TJ this spectrum 
is dominated by strong incommensurate wings below the 
resonance which are in fact predicted to be much more 
pronounced than in the case of the p-type supercon- 
ductors. This is precisely opposite to the experimental 
findings where the incommensurate fluctuations are pro- 
nounced in the p-type systems, but completely absent in 
the n-type superconductor! 

Including a small q-modulation of the form U{q) = 
Uo — AU[cos{qx) + cos{qy)] as used recently^^ with a rel- 
ative modulation AU/Uq = 0.1 to obtain a slightly bet- 
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FIG. 6: Comparison of the resonance in absolute units with 
FL/RPA calculations for optimally hole-doped YBCO and 
electron-doped PLCCO. (a) Local susceptibility in absolute 
units for optimally doped YBCO at 10 K from Ref. Q. The 
solid blue line is the calculation based on RPA model scaled to 
match the experimental data, (b) Local susceptibility in both 
the normal and superconducting state for PLCCO obtained 
from converting the raw data of Ref. [l^ to absolute units. 
Solid lines are guide to the eyes. The dashed lines represent 
the results of the FL/RPA calculations with the same scale 
factor as used for YBCO. Note that the theoretical values are 
about 6 times smaller than the experimental results. 



ter quantitative agreement with continuously improving 
neutron scattering data on optimally doped YBCO, does 
not lead to significant improvements but only to a small 
change of the incommensurability of the wings. 

The only way to repair this gross inconsistency is 
by invoking a U{q) which sharply peaks at qaf- Re- 
cently, it was argued^^ that by taking a full momen- 
tum dependent four-point vertex (C/q 0), U{q} = 
—J[cos(qx) + cos((7j,)]/2, the incommensurate wings can 
be suppressed. This strongly momentum-dependent form 
of the four-point vertex peaking at qaf was motivated 
by the proximity of the superconducting and (commen- 
surate) antiferromagnetic phases. However, such a form 
of the four-point vertex is clearly unphysical since it cor- 
responds to a nearest neighbor exchange whereas the on- 
site Coulomb repulsion which is known to control the 
Mottness in the copper oxide planes is completely ig- 
nored. In contrast to the i J-model, in the quasiparticle 
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FIG. 7: Comparison of the magnetic excitation spectrum 
x"{ci,Lj) along the [H,H] direction in the vicinity of qap = 
(1/2, l/2)r.l.u. resulting from the FL/RPA calculations with 
neutron scattering data^^ on optimally doped PLCCO {Tc = 
24K) measured at T = 2K. The very strong incommensurate 
wings predicted by the calculations highlight the failure of the 
FL/RPA approach. 



picture used here double occupancies are not projected 
out. Moreover, using our set of quasiparticle parameters 
(Tab. IJ), this would imply an effective superexchange of 
J = 500 meV which is about 5 times bigger than in the 
parent undoped compounds^l. While this is obviously 
unphysical, the value J — 854 meV taken in Ref. 26 is 
even much bigger. 

The reason why incommensurate wings at low ener- 
gies appear generically within the FL/RPA approach for 
any realistic set of parameters both for p-type and n- 
type materials is actually a very generic one, rooted in 
the assumption that there is a direct relation between 
the free particle-hole and the magnetic spectrum. Within 
this framework, the RPA response x"{(l, oj) for any realis- 
tic form of J7(q) basically reflects the momentum depen- 
dence of the gap of the ph-continuum nicely seen in the 
bare Lindhard function XQ{q,uj) (Fig. [5]). The supercon- 
ducting d-wave gap is close to its maximum for particle- 
hole pairs separated by qaf and goes continuously down 
if we move away from the antiferromagnetic to incom- 
mensurate wavevectors separations (see Figs. Ill9p . The 
gap of the ph-contiuum closes at the incommensurate 
wavevectors connecting points of the Fermi surface coin- 
ciding with the nodes of the c?-wave gap. 

Although the momentum dependence of the gap of the 
particle-hole continuum looks very similar in the p-type 
and n-type case, a crucial difference becomes apparent 
when comparing the distribution of the spectral weight 
Xg(q, w). Whereas in p-type YBCO spectral weight is 
accumulated at qaf the intensity in the close vicinity of 
the antiferromagnetic wavevector is strongly suppressed 
in n-type PLCCO. On the other hand the spectral weight 
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FIG. 8: Bare Lindhard functions Xo(q,'^) of YBCO and 
PLCCO in the superconducting phases calculated with band- 
structure and d-wave gap parameters listed in Table U] 
Whereas the momentum dependence of the gap of the ph- 
continuum looks very similar for the p-type and n-type ma- 
terial, the distribution of spectral weight is completely differ- 
ent. For the p-type, spectral weight is accumulated at qAp 
whereas for the n-type a lot of intensity has shifted from qaf 
to incommensurate wavevectors. 



at incommensurate momenta is strongly enhanced in the 
n-type compound (Fig. [5]). 

The reason for the reshuffling of the weight in the bare 
Lindhard function Xq (q, w) from Qaf to incommensurate 
wave vectors in going from p- to n-type superconductors 
is simply related to the number of particle-hole pairs con- 
tributing to the magnetic response. In Fig. [2]we compare 
the normal state dispersions, Fermi surfaces, and Fermi 
velocities of YBCO and PLCCO. Whereas in YBCO the 
saddle-points in the band structure responsible for the 
van Hove singularities at the antinodal points are very 
close to points on the Fermi surface separated by Qaf, in 
PLCCO the bands are very steep at points connected by 
qaf and a nesting of the Fermi surface for incommensu- 
rate wavevectors in regions of very flat bands give rise to 
the drastic spectral weight enhancement of the wings. 



V. DISCUSSION AND CONCLUSION 

To summarize, by combining experimental measure- 
ments of the quasiparticle and dynamical magnetic prop- 
erties we have demonstrated that in a n-type cuprate 
superconductor the magnetic excitations to be expected 
from a weakly interacting Fermi-gas are inconsistent with 
experimental observations. In particular, we have per- 
formed tunneling experiments on the same sample of 
PLCCO showing a magnetic resonance in the supercon- 
ducting phase^^ and used ARPES datgfci^ on a similar 
sample to extract normal-state band structure and d- 
wave gap parameters. The neutron scattering raw data 
on the magnetic resonance^^ we have converted to abso- 
lute units by normalization to acoustic phonons. 

Using the detailed information on the quasiparticles 
obtained from the ARPES and tunneling experiments we 
have calculated the expected magnetic excitation spec- 
trum within the conventional FL /RPA framework which 
assumes that there is a direct relation between the free 




FIG. 9: Comparison of the normal state dispersions, Fermi 
surfaces (left panel) and Fermi velocities (right panel) of 
YBCO and PLCCO. 



particle-hole and the magnetic spectrum. The compar- 
ison of the theoretical results with the magnetic fluc- 
tuation measured in inelastic neutron scattering shows 
that the fermiology approach fails to explain the mag- 
netic fluctuations. 

Since the magnetic resonance of PLCCO is located at 
an energy near the gap of the particle-hole continuum as 
confirmed by our tunneling experiment, it is difficult to 
explain it as a bound state within the FL /RPA approach. 
Consequently, within the FL/RPA framework we obtain 
an almost step-like feature rather than a symmetric res- 
onance peak seen in experiment and underestimate the 
spectral weight of the resonance by an order of magni- 
tude. Additionally, taking the temperature dependence 
of the gap measured by our tunneling experiment, the 
FL/RPA approach predicts a very strong temperature 
dependence of the resonance well below Tc inconsistent 
with the experimental observation. The failure of the 
fermiology framework is highlighted by the incommen- 
surate wings which from the theoretical calculations are 
expected to be much more pronounced than in the p-type 
case whereas they have never been observed in electron- 
doped superconductors. 

Within the FL/RPA approach such incommensurate 
wings in the magnetic response of a d-wave superconduc- 
tor appear generically for any physically reasonable set 
of parameter, both in the p-type and n-type case. This 
finding is quite robust and does not depend on details of 
the bandstructure. However, the different forms of the 
quasiparticle dispersion and of the d-wave gap of PLCCO 
compared to YBCO gives rise to an additional reshuf- 
fling of spectral weight in the free particle-hole spectrum 
from the antiferromagnetic to incommensurate wavevec- 
tors leading to an enhancement of the wings and an ad- 
ditional intensity loss in the close vicinity of Qaf- 
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The drastic failure of the fermiology approach for the 
n-type case opens the question whether the apparent 
agreement for the p-type superconductors is just coin- 
cidental. Since incommcnsureate wings are generically 
expected within the FL/RPA approach it is not surpris- 
ing that one finds a reasonable agreement up to the reso- 
nance energy. However, this approach cannot explain the 
upper branches of the hourglass spectrum seen in vari- 
ous experiments. More severely, the fermiology interpre- 
tation can neither account for the anomalous properties 
of the normal state which is known to be a non-Fermi 
liquid nor for the persistence of the resonance and the 
hourglass above Tc in the underdoped regime. 

On the other hand, above the spin gap the mag- 
netic excitation spectra of superconducting YBCO^i^ and 
La2-2;Sr2;Cu04^ are remarkably similar-^'^* to that found 
in stripe ordered Lai.875Bao.i25Cu04i suggesting that 
the magnetic fluctuations in the p-type superconductors 
correspond to fluctuating stripes competing with super- 
conductivity. Theoretically, the hourglass spectrum char- 
acteristic for both, stripe ordered and superconducting 
p-type cuprates, has been obtained in various models 
for static stripes^^ but also in a phenomenological lat- 
tice model for thermally fluctuating, short-ranged stripe 
order—. 

Whereas the magnetic fluctuations in the p-type 
cuprates seem to reflect the competition between su- 
perconductivity and incommensurate, Mott-like antifer- 
romagnetism (the "stripes"), incommensurate fluctua- 
tions have never been observed in n-type superconduc- 
tors suggesting instead a competition with commensurate 
antiferromagnetisroii. 

To conclude, by combining experimental measure- 



ments of the quasiparticle and dynamical magnetic prop- 
erties we have demonstrated that in the n-type cuprate 
superconductor PLCCO there is no relation whatsoever 
between the magnetic excitations to be expected from 
a weakly interacting Fermi-gas and the magnetic fluc- 
tuations observed experimentally. This demonstrates 
that the magnetic fluctuations actually correspond with 
highly collective motions which likely reflect the quan- 
tum competition between superconductivity and strongly 
coupled antiferromagnetism. The challenge for the the- 
orist is to explain how this system manages to simul- 
taneously support conventional looking fermionic quasi- 
particle excitations and highly collective order parameter 
fluctuations. 
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